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Using the standard concepts of free random variables, 
we show that for a large class of nonhermitean random ma- 
trix models, the support of the eigenvalue distribution follows 
from their hermitean analogs using a conformal transforma- 
tion. We also extend the concepts of free random variables 
to the class of nonhermitean matrices, and apply them to 
the models discussed by Ginibre-Girko (elliptic ensemble) and 
Mahaux-Weidenmiiller (chaotic resonance scattering). 



1. Introduction 

The distribution of eigenvalues of complex and large 
random matrices is of relevance to a variety of physical 
problems. Non-hermitean random matrices appear natu- 
rally in the evolution of dissipative quantum many-body 
systems Jlj], in quantum chaotic scattering ||, in quan- 
tum optics fl , and possibly in quantum chromodynamics 

The distribution of complex energies of unstable quan- 
tum systems in the framework of random matrix mod- 
els, has been investigated in some details by Sokolov and 
Zelevinsky ||, following on the original work of Wei- 
denmiiller and coworkers Alternative and extended 
analyses can be also found inH using the replica method, 
and in |^| using the supersymmetric method. For large 
matrices, a structural change in the eigenvalue distribu- 
tion within the complex plane was observed in the case of 
strong nonhermiticity. The level density has been used to 
assess the statistical distribution of the resonance widths 
in chaotic scattering with large scattering channels |7j, 
and the time delay in chaotic scattering ||]. The lat- 
ter may be of relevance for a quantitative assessment of 
quantum chaos 

In this paper, we would like to show how to evaluate in 
a straightforward way the supports for the level density 
as well as the eigenvalues distribution for a large class of 
nonhermitean and random matrices using the concepts 
of free random variables as developed by Voiculescu flCj l 
and popularized by Zee [jll|. In section 2, we outline 
the general definitions for the standard S and R trans- 
formations, as known for hermitean random matrices. In 



section 3, we discuss the resonance matrix model as an 
example of a nonhermitean random matrix model. In 
section 4, using the concepts of S and R transforms we 
derive the spectral density of its hermitean analogue. In 
section 5, we show how the support for a nonhermitean 
spectral density follows from the hermitean one using a 
conformal transformation. Our arguments are applied 
to Ginibre-Girko's (l2) and Mahaux-Weidcnmullcr's p"3| ] 
random matrix models. In section 6, we generalize the 
concept of R-transforms to the nonhermitean case, and 
use it to analyze the resolvent of nonhermitean random 
matrix models. Our conclusions are summarized in sec- 
tion 7. 

2. R and S Transformations 

Addition and multiplication of free random variables 
can be assessed using R and S transformations p0[ . 
Specifically, the R-transformation is additive -R1+2 = 
i?i + i?2, and the the S-transformation is multiplicative 

For the sum, if we were to define a Blue's function 
B(z) jllj as the functional inverse of a Green's func- 
tion G(z), that is B(G(z)) = G(B(z)) = z, then the 
R-function is simply R(z) = B(z) + 1/z. Physically, 
the R-transform is some pertinent self-energy in the pla- 
nar approximation The additive property of the 
R-transform implies the addition law 1 10[ Tl| ] , 



B 1+2 (z) =B 1 {z) + B 2 {z) 
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(1) 



for the Blue's functions. Hence, the problem of finding 
the spectral distribution of the sum of two free random 
matrices is linear, and follows from the simple algorithm: 
Find Gi and G2, construct their functional inverses B\ 
and B2, add them through ([!]), and invert S1+2 to get 
Gi+2- The spectral density of the sum is the discontinu- 
ity of G1+2 along the real axis. The zeroes of B'{z) = 
characterize the end-points of the spectral density , 
and reflect structural changes in the underlying spec- 
trum. 

For the product, if we were to define x{ z ) as a solution 
of @ 
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then the S-function is simply 
S(z) = 



-\ = z 



(2) 



(3) 



The S-transform of the product of two free random ma- 
trices is the product of their S-transforms that is 6*1*2 = 
5*1 5*2. Given 5*1*2, the resolvent Gi*2 follows through (j|) 
and (H) in reverse order. 

3. Resonance Scattering Model 

To make part of our subsequent discussions clear, we 
will use the nonhermitean random matrix model intro- 
duced by Mahaux and Weidenmuller for resonance scat- 
tering as inspired by Weisskopf and Wigner effective 
hamiltonian |l4|], to illustrate some of our assertions. In 
brief, a quantum system composed of (N — M) closed 
and M open channels can be described by the effective 
scattering Hamiltonian 



H = H E 



igT 



T = AA T 



(4) 



which is N x N dimensional. A is an asymmetric N x 
M random matrix, and g an overall coupling. Unitarity 
enforces the form of T used in (||). Wigner's condition 
implies g < 0. The matrix elements A\ characterize the 
transition between the (N—M) -internal channels and M- 
external channels, and are assumed to be independent of 
the scattering energy. For M = 0, the Hamiltonian is real 
and the spectrum is bound. For M^O, the Hamiltonian 
is complex, with all states acquiring a width. 

The resolvent G(z) associated to (0) is defined as 



(5) 



with complex valued poles, reflecting the energy and 
width of the resonance states. Its imaginary part is re- 
lated to the total life-time (time-delay) in the scattering 
process. The averaging in (||) corresponds to the GOE 
ensemble for Hp, with the transition matrix elements Af. 
treated as independent Gaussian variables . The sim- 
pler case ("f-case" Hfl)), where the matrix elements are 
fixed, will be briefly discussed later. In its domain of an- 
alyticity D, G{z) is a holomorphic function of z. In the 
complement D, G{z) is in general nonholomorphic, with 
a nonvanishing spectral distribution [flEf 



IdG(z) 
v(z,z) = — 

IT OZ 



(6) 



The latter provides for a statistical analysis of the char- 
acteristics of the resonances in chaotic scattering for a 
large number of channels J?]||[l(§ • It can also be used to 
discuss Dicke superradiance in quantum optics ^ . 

4. Hermitean Case 



First consider the case of an N x N real symmetric 
product matrix = AtAf. In the large ./V limit, the 
spectral function is connected to the one of the random 
NxN matrix A as 



NxN 



(z = A 2 ) 



IAI 2tt Jl 



(7) 



The numerator of (Q) is just Wigner's semicircular 
law p7[ ]. The case of rectangular N X M matrices with 
TV, M — > 00 but m = M/N fixed, follows by truncation 
using the projector pl| 



P = diag(l 1 _^ 1 l,0 : __0) 

M N-M 



(8) 



that is T = V s P. We recognize immediately the prob- 
lem of "multiplication" of the random matrix V s by the 
deterministic projector P. Following the multiplication 
algorithm we construct below the resolvent for 

the product. First, we construct the resolvent for the 
projector 



Gp(z) = — Tr— ^~ 
V ; N z-P 

The S-transform is 

Sp(z) -- 



1 , A 

= m h (1 — to) — 



1 + z 



TO + Z 

The resolvent for the square of the Gaussian is 



G A 2 = Ul - w i - - 



(9) 



(10) 



(11) 



Its discontinuity along the real axis reproduces (f7j). The 
corresponding \ and S transformations are 



Xa 2 



(1 + z) 2 ' 
so the product matrix has 



S A 2 = 



l + z 



S Syj2 S P 



1 

TO + Z 



and 



Xa 2 *p 



(1 + z){m + z) 



Inverting (14), and inserting z(x) m to 
resolvent 



G{z) 



1 — TO 1 

2z + 2 




(12) 

(13) 
(14) 

we get for the 
(15) 



The spectral density follows from the discontinuity of 
( |T5| ) along the real axis in the form 
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(16) 



This result was first obtained in [fL9[ , using other meth- 
ods. The first term, originating from the pole 1/z repre- 
sents (N — M) zero modes of the matrix V. 

Now we could "add" the random hamiltonian Hr. It 
is convenient first to solve the hermitean problem Hh — 
Hu+gT. This problem is immediately solved using the R 
transformation (Blue's function), since the resolvent for 
GOE is known to be 



G{z) = Vz 2 -4) 



(17) 



and the resolvent for T is given by (|15j) . The correspond- 
ing Blue's functions are straightforward to find 



Bi(z) 
B 2 (z) 



1 

z + - 

z 

mg 
1 - gz 



(18) 



where we have reinstated the coupling g. The labels 1(2) 
refer to iJ^(r)) respectively. From the addition law fol- 
lows the resolvent G of the sum Hh, as a solution to a 
cubic equation (Cardano class) 



1 



mg 



G(z) l-gG(z) 



G{z) 



(19) 



Out of the three solutions to this equation, we choose 
uniquely the one respecting the positivity and normaliz- 
ability of the spectral density. Structural changes in the 
spectral density can be easily read out from the the zeros 
of B'(z) = (end-points), where B is the Blue's function 
for the sum. This equation is fourth order (Ferrari class) 



1/z 2 



mg* 



(1-gz)- 



+ 1 = 



(20) 



and defines all the endpoints of the spectrum Ai =B(zi), 
(i = 1, 2, 3, 4). A structural change (usually a phase tran- 
sition in the underlying system) happens when two arcs 
[Ai, A2] and [A3, A4] merge into the one (i.e. A2 — A^)^ 
The condition for the structural change (merging point) 
follows from the zero of the discriminant of the quartic 
equation (po|), 

(1 - (1 - m)g 2 ) 3 + 27g 4 m = , 



with the solution 

9* (m) 



1 



(1 - ^m) 3 



(21) 



1 In the present case, the spectrum is not even in contrast to 
the chiral cases considered in [ECl. 



Before closing this section, let us mention that the sim 
pier f-case [[ 
strained by 



where the matrix elements A k are con- 



J2 A a k A b k =S ab (a,6=l,...,M) 



(22) 



k=l 



instead of being Gaussian, can be analyzed using similar 
methods. In a reaction process, the constraint (|2^) ex- 
cludes direct reactions. In this case, the eigenvalues of 
the matrix gT are either zeros (since A is rectangular) 
or equal to g (because of (22)). The resolvent for T in 
this case has the form of a projector, with a fraction m of 
eigenvalues equal to g, and the remaining fraction (1 — m) 
equal to zero 



Gt(z) = 



1 — m 



(23) 



The addition law for the Blue's functions yields the cubic 
Pastur's equation Ell] 



G 



1 



/ 



G 



f 



-G 



f-9 



(24) 



Repeating the analysis for the end-points of the previous 
example shows that a structural change in the spectrum 
takes place at 



(25) 



gl f (m) = [s^l -m + ^m]' 



While the first order terms in ym in (|21|) and ( J25 
equal the remainder is not. Note that ( pi"] ) diverges 
as m — » 1, while (^5|) is symmetric around m = 0.5, with 
gl varying between 1 and 4. 

5. Conformal Mapping 

Under the substitution g — > ig the matrix ensemble 
becomes complex, taking us to the resonance scattering 
model. As a result, the eigenvalue distribution is valued 
in the complex plane. For large z, the resolvent G(z) is a 
holomorphic function of z, whose form follows from jl9| ) 
through g — > ig. As z is decreased, G(z) will in general 
blow up (zero denominator), a signal that the function is 
no longer holomorphic. 

The boundary between the holomorphic and nonholo- 
morphic solutions in the z-plane, can be derived very gen- 
erally using a conformal transformation that maps the 
cuts of the hermitean ensemble onto the boundary of its 
nonhermitean analogue. Indeed, consider the case where 
a Gaussian random and hermitean matrix R is added to 
an arbitrary matrix M. The corresponding Blue's func- 
tion is 

B R+M (u) = B R {u) + B M {u) --= B M (u) + u (26) 

u 

Substituting u — > Gm{z) we get 
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B 



R+M 



(G M (z)) = z + G m {z) 



(27) 



Let wbea point in the complex plane for which Gm(z) = 
Gr+m{w)- Using the definition of Blue's functions, it 
follows that w is located at 



w 



z + G M (z) 



Now, if we were to note that in the holomorphic domain, 
the Blue's function for the Gaussian and nonhermitean 
ensemble is Bm = —z + 1/z, then 



B 



R+iR 



= B f 



D 1 _ 1 

t>m — - 



(29) 



The antihermitean Gaussian nullifies the contribution of 
the hermitean Gaussian (the R-function in this case is 
zero). This observation together with (p8| ) allows for re- 
lating a hermitean resolvent to a nonhermitean one in the 
holomorphic domain. Indeed, using (^8|) and the present 
observation, it follows that 

w = z- 2G R+M (z) G iR+M (w) = G R+M (z) (30) 

throughout the holomorphic region. Similarly, 

w = %{z - 2G R+M {z)) G R+iM {w) = -iG R+M (z) (31) 

after a simple rotation of the real and imaginary axes. 
The transformation ( |3l"| ) maps the cuts of the hermitean 
ensemble onto the boundaries between the holomorphic 
and nonholomorphic regions of the nonhermitean ensem- 
ble, as we now illustrate. 

• Elliptic Distribution 
Consider first the simple case where H = H R + gH R , 
with H R random Gaussian. Since the Blue's function in 
the holomorphic region is B = 1/z + {1+ g 2 )z, then the 
resolvent (inverse) is just 



9 2 G, 



The two solutions to 
1 



t) are 



G(z) = 



2(1 



±^z 2 -A(l + g 2 ) 



(32) 



(33) 



( |33| ) are holomorphic everywhere in the z-plane, except 
on the cut C = [—2^Jl+g 2 , 2^/l+g 2 } along the real axis. 
C is the support for the spectral density associated to H 
(Wigner's semi-circle). 

In the holomorphic region, the non-hermitean case fol- 
lows by setting g to ij yielding B = 1/z + (1 
Hence 



l 2 )z- 



w 



1 



7 2 G a . 



(34) 



for the resolvent. Using the hermitean solution ( p3| ) and 
the mapping (BCj), we can map the hermitean cut C onto 



the boundary delimiting the holomorphic region for (|34|), 
that is 



w 



— {g 2 z ± V* 2 -4(l + fl »)} 



(35) 



(28) with z — t ± iO and t in C. Equation ( |35| ) span an ellipsis 
with axes 2/^/l+g 2 and 2g 2 /y / l+g 2 . For g 2 = 1 this is 
just Ginibre's circle^]. 





g=-0.1, m=0.25 











g=-l,m=0.25 
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FIG. 1. The evolution of the boundary (thin solid line) 
for various couplings at fixed ratio m=0.25. The dots are 
the numerical eigenvalues generated from an ensemble of 50 
matrices that are 100 by 100. The solid bars are the positions 
of the cuts following from (|20|) with g — > ig. 



• Resonance Scattering Model 
In the same spirit, the cuts associated to H R + gT, and 
given by (20), can be mapped onto the boundary dD 
delimiting the holomorphic region for the nonhermitean 
ensemble H R + igT. The mapping ( ^l| ) requires the so- 
lution to ( |19|) for the hermitean ensemble. Indeed, after 
rewriting ( |19| ) in terms of w and z, explicitly inserting the 
form w — x + iy and choosing z = t (where real t scans 
the cuts of (|l^)) one arrives at two coupled equations for 
x, y and t. Eliminating t, the equation of the boundary 
reads 



4m 

gy 



gy 



m 

y 



(36) 



2 Here for GUE the scale is set by (3 — 2 yielding the radius 
of the disc to be V2. In the case of GOE, where the implicit 
scale is /3 = 1, the radius is 1. 



4 



in agreement with pf. 

The results are displayed in Figure [l] for different cou- 
plings g and a fixed ratio m — 0.25. The locations of the 
cuts following from the holomorphic resolvent through 
© with g — > are found to be strongly correlated 
with the shape of the envelopes following from the non- 
holomorphic resolvent. The envelopes can be regarded 
qualitatively as a smearing of the cuts. For small cou- 
plings, say \g\ — 0.1, the vertical cut lies far on the x — 
axis in the positive y half-plane. It does not show up 
in our first figure. For larger couplings, say \g\ = 1, it 
emerges from the half-plane, and gives rise to a bulge for 
\g\ = 4. Above the critical value \g\ = 4.44 following 
from (pl|), the vertical cut drives the upper island up, to- 
wards y — — oo for \g\ — > oo. At very large couplings, say 
\g\ = 100, the islands are indistinguishable from the cuts. 
The stability of the lower island composed essentially of 
long-lived states (small width states) is ensured by the 
zero modes of T, as originally suggested by ||||. 

Repeating the calculation above for the f-case, using 
the same transformation (pTJ) with a resolvent G now 
taken from the solution of (|24|), yield the following bound- 
ary 



4m 

gy 



y 



y 
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(37) 



A diagrammatic proof of this relation will be given else- 
where |p2f . The present method can be used to analyze 
the Green's functions of hermitean and nonhermitean 
random matrix models in the entire z-plane. We now 
illustrate this assertion for the two cases discussed above. 

• Circular Ensemble 
Consider once more the case Hr + iHr. The generalized 
Green's function for the hermitean Gaussian ensemble 
Hr follows from a straightforward generalization of Pas- 
tur's equation [pTJ (resummation of the rainbow graphs) 



1 



Z-T, 



through the substitution z — > Z. Hence, 

G + i = Z 
G 

The corresponding Blue's function is 
B R (A) = j+A 



(41) 



(42) 



(43) 



which is to be compared with (|18j) . For the nonhermitean 
Gaussian ensemble, the new version of Pastur's equation 
is 



between the holomorphic and nonholomorphic solutions. 

We note that a similar construction can be used to 
analyze the nonhermitean chiral random matrix model 

i- 

Let us finally mention, that the technique described 
here does not allow directly to go inside the islands of 
non-analyticity, where the Green's functions are nonholo- 
morphic. For that we need to generalize the concepts of 
free random variables for nonhermitean random matrices 
as we now discuss. 



6. Blue's Functions Revisited 

To discuss the generic case of nonhermitean random 
matrix models, we need to generalize the concept of R 
transformation. For that, consider the addition law (Q), 
rewritten in the equivalent way (z — > G\ + 2{z) = G), 



z = Bi(G) + B 2 (G) - 



G 



(38) 



The generalization to nonhermitean random matrices 
amounts to defining 2x2 Green's functions G, and Blue's 
functions B, such that 



B(G) = Z = 



z 
z 



The addition law (R-transform) becomes 
Z = B 1 {G) + B 2 (G) - 1 



(39) 



(40) 



Go 



-1 1 
1 -1 



= E 



(44) 



where the multiplication is meant componentwise. The 
extra matrix sign follows from a sign flip in the "gluon- 
propagator" (nonhermitean case), while summing over 
the rainbow graphs |22]]. The resulting Blue's function is 



«^)-3 + (i:Hi-! 



(45) 



The Green's function for the sum comes from (f40|). In- 
deed, if we were to define 



G = 



a b 



then © together with © and © give 



a b 
b c 



(46) 



—a b 
b -c 



The equation for the off-diagonal element 







det£ 



2b 



(47) 



(48) 



has two types of solutions. The case b = corresponds to 
the holomorphic solution a = G = 1/z discussed above, 
since a in G is just the resolvent (|5|). The nonzero case 



5 



with detQ = 1/2, gives the nonholomorphic (here anti- 
holomorphic) solution a = G = z/2 in the complemen- 
tary part of the complex plane, with a uniform eigenvalue 
distribution u(z, z) — 1/2tt as expected from (pi) . This is 
just the result derived by Ginibre and Girko [12]. 

• Random Scattering Model 
Here we consider, for convenience, the ensemble of ma- 
trices 



H 



-gAA 



(49) 



with again A an N x M random complex matrix, Hr a 
hermitcan random matrix. The hamiltonian H is related 
to hamiltonian (||) by multiplicative factor i, resulting in 
rotating the axis of the eigenvalue plane by n/2. The 
generalized Blue's function for —gAA T follows from JI§| ) 
(second equation) through the substitution z — > A, g — > 
—g. Hence, 



B. 



gAA 



■(A) = - 



mg 



1 



1+gA A 
Using the addition formula we get 



z 
z 



1 + gG 



i 
G 



-1 
1 



1 

-1 



(50) 



(51) 



This equation admits a holomorphic and nonholomor- 
phic solution. The former is in agreement with the one 
discussed above, while the latter reads (z = x + iy) 



G 



vy 
2 



rn 

x 



9 



1 - gx 



(52) 



in agreement with |7j . The distribution of eigenvalues has 
a support only in the nonholomorphic region and follows 
from (^). Here, the boundary between the holomorphic 
and nonholomorphic solution follows from the vanishing 
of the off-diagonal elements of Q in (^), reproducing 
exactly (|36|). 

Deriving the resolvent for the f-case is straightforward, 
since the problem reduces to finding the solution to the 
generalized Pastur equation J22| 



1 



with 



Gf = 



Z-Gf 



-1 
1 



Of 



Gf 



(53) 



(54) 



The nonholomorphic solution for the resolvent (|^) reads 
1 ( m — 1 



G i 



W 
2 



g + x 



rn 

x 



2x-g 



(55) 



in agreement with 



7. Conclusions 



We have shown how the concepts of addition and mul- 
tiplication of free random variables could be used to an- 
alyze some problems related to the hermitean and non- 
hermitean random matrix models, in the large N limit 
and for the case of strong nonhermiticity. For the ran- 
dom scattering model, our level and width distributions 
are in agreement with those discussed in Q for strongly 
overlapping resonances. In this sense, we have not much 
to add to their general physical discussion. 

Our approach extends the work of Voiculescu and Zee, 
and to our knowledge is new. In particular, we have 
shown that the supports of complex eigenvalues for the 
nonhcrmitean ensemble follow from a simple conformal 
transformation on the cuts of its hermitean analog. Our 
approach offers a simple alternative to the replica and 
supersymmetric methods for a variety of random matrix 
models in the large N limit. A comprehensive compar- 
ison, including the issue of 1/N corrections, is unfortu- 
nately beyond the scope of this work. 

The resonance scattering model discussed here and the 
solutions presented above, could be used to discuss some 
issues related to chaotic behavior in complex systems 
such as the instanton liquid model, conductance in quan- 
tum dots, microwave cavities or nuclear dynamics. We 
hope to address some of these issues next. 
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